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Abstract 

First, we show that, if there are no bound states, we can express the q.m. 
third cluster - involving 3 and fewer particles in Statistical Mechanics - as a 
formula involving adiabatic eigenphase shifts. This is for Boltzmann statistics. 

From this q.m. formulation, in the case of purely repulsive forces, we 
recover, as h goes to 0, the classical expressions for the cluster. 

We then discuss difficulties which arise in the presence of 2-body bound 
states and present a tentative formula involving eigenphase shifts and the 2 
and 3 body bound state energies. We emphasize that important difficulties 
have not been resolved. 



Statistical Mechanics 

In equilibrium Statistical Mechanics ALL wisdom derives from the partition func- 
tion! Here, we need the logarithm of the Grand Partition function Q: 



In Q = z Tr{e 
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which when divided by V, gives coefficients which are independent of the volume, 
when the latter becomes large; we call them bi. The fugacity z equals exp{^/ kT), 
where fi is the Gibbs function per particle, k is Boltzmann's constant and T is the 
temperature; f3 = 1/kT. We can then write for the pressure and the density 

p/kT = {l/V) \nQ = Y.hiz' 

I 

N/V = p = Y.lbi z' 
I 

The fugacity can then be eliminated to give the pressure in terms of the density. 

p/kT = p-\ 

The coefficients of the second and higher powers are called the virial coefficients. 
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Crucial Step 

For this work we extract the Boltzmann part of the traces: we write 
Tr{e~^"^) = —Trace'^ie-^^") + Exchange Terms 

We can then write for the Boltzmann 63: 

63 = (3!y)-iTrace^[(e-'^^3 - e'"^-') - 3 (g-'^^^^+^i) - e''^^^)] 

where I have made use of the Boltzmann statistics to express the answer in terms 
of 3-body traces. 



Adiabatic Preliminaries 

For the 3 particles of equal masses, in three dimensions, we first introduce center of 
mass and Jacobi coordinates. We define 

V- - r2) e = (^)'^'(^^^ - ^1^) I + + ra) 

where, of course, the n give us the locations of the 3 particles. This is a canonical 
transformation and insures that in the kinetic energy there are no cross terms. 

The variables ^ and 77 are involved separately in the Laplacians and we may 
consider them as acting in different spaces. We introduce a higher dimensional 

f 

vector p={_^) and express it in a hyperspherical coordinate system (p and the set 

of angles O). If we factor a term of p^/^ from the solution of the relative Schrodinger 
equation, i.e. we let t/j — (f)/ p^^'^, we are lead to: 



92 2mE 



0(p,Q) = O 



where 



and m is the mass of each particle, E is the relative energy in the center of mass. 
Vfj is the purely angular part of the Laplacian. We now introduce the adiabatic 
basis, which consists of the eigenfunctions of part of the Hamiltonian: the angular 
part of the kinetic energy and the potential. 

H,Biip,n)^Aeip)Bi{p,n), 

where i enumerates the solutions. 



2 



Using this adiabatic basis, we can now rewrite the Schrodinger equation as a 
system of coupled ordinary differential equations. We write 

and obtain the set of coupled equations 

+ = 0, 

where k'^ is the relative energy multiplied by 2m/ and we defined: 

Di^M = ldnB*{n,p)—Be>{n,p). 

We note that 

i)«. = ^(C„,) + (C=),,, 

The Phase Shift Formula 

When there are no bound states, we may write 



Tr^{e-^"^)^ jdp jdkY,'4^\k,p){'il)\k,p)y e 

where we have introduced a complete set of continuum eigenfunctions. Expanding 
in the adiabatic basis, we obtain 

rr^(e-^^^) = jdp jdkY^mKp)? e-^^&^"\ 

where we note that we have integrated over the angles and taken advantage of the 
orthogonality of our i?/'s. Wc integrate from to oo. 

We now return to our expression for 63 and proceed as above, but drop the tildas, 
to obtain: 



3^^ 
2 



Jdke-^^'^ jdpY,m\' - \<i>i,\') - - \4>i,\')i 



where we have evaluated the trace corresponding to the center of mass. The ampli- 
tudes (p\ correspond to ifa, (p[ to H2 + Ti and amplitudes with a zero belong to the 
free particles. The thermal wavelength At is defined as h/ \/2wrni^. 
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We now make use of a trick to evaluate the p integrals. We first write 

/ mk, p)\' dp = lim / <t>\{K p)m', p) dp 



and then, and there is the trick, 







j^im^p) m\p))dp^ 



evaluated at p — pr, 



I.e. our identity is: 



X ^ d 

'^'dp 



dp dp 



+ {k'-{k'f)Y.Mk)Mk') 

i 

e,e' ^P 

and we integrate with respect to p. Using then the fact that (p goes to zero, as p 
itself goes to zero, and that C decreases fast enough for p large, we are left with the 
expression displayed earlier (that of our 'trick'). 



We now put in the asymptotic form of our solutions, oscillatory solutions valid for 
Pmax large, and use I'Hospital's rule to take the limit as k' k. 
The solutions are: 

(f)} {kpy/'^Ce,i [cosSi JK+2{kp) - smSi NK+2{kp)] 

where the order K is one of the quantities specified by £. Inserting this into our 
integrals we find that 

rPmax • 2 Id- 1 

/ \<Pe{k)\ dp -—5\k) + -pmax + osc. terms 
^ Jo n dk TT 

and, thus, that 

mk)\' - |0^,o(^)r) dp ^ --^^^k) + osc. terms 

We let Pmax go to infinity, and the oscillating terms - of the form sin{2kpmax + ■■■)" 
will not contribute to the subsequent integration over k. A partial integration now 
gives us our basic formula. 
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where 

i 

The first S arises from comparing three interacting particles with three free particles. 
The second 5 arises when a 3-body system, where only two particles are interacting 
(one particle being a spectator), is compared to three free particles. 



Classical Limit 

The idea behind our WKB treatment of our equations, is to argue that when the 
potentials change slowly - within oscillations of the solutions - then the adiabatic 
eigenfunctions will also change slowly and we can neglect their derivatives. Thus 
we will obtain uncoupled equations with effective potentials (the eigcnpotentials 
A^(p)). We then proceed with these in a more or less conventional WKB fashion. 
Let us assume, here, one turning point po- 

The phases can now be obtained by considering simplified forms of the asymp- 
totic solutions for the Let us denote them as (f)„. The phases will then be 



5,r.(K + 2)--kpo + [Jk^-A.-^-k]dp 

Inserting our expression for S,, into dk k S^{k) exp{— X^k^ / Att) and interchanging 
the order of integration (p and k) we obtain: 

Summing now over c, we can rewrite the exponentials as traces: 



Tra.e«{exp[-| (A(p) + ±)] - exp[-| ^]} 



where A is the operator (matrix) which yields the diagonal elements A^, and /C^ 
the operator which yields the eigenvalue when the interaction is turned off (and 
therefore takes on the diagonal values (i^+2)^— |, associated with the hyperspherical 
harmonic of order K). The trace is restricted so as not to involve p. 
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In another key step, we switch to a hyper spherical basis. We note that A is related 
to (2m/^^) V + K? / by a similarity transformation and an orthogonal matrix U. 
Substituting in the trace, we lose the U and obtain 

X2 I 1 \2 r-2 I 1 

^"l-P(-/3l'-^^)-exp(-^^)l 

We write the exponential as a product of 2 exponentials, disregarding higher order 
terms in %. Introducing eigenkets and eigenbras which depend on the hyperspherical 
angles, we write the trace as: 

r \2 1^2 I 1 

dQ<n\ exp(--^ > {exp[-pV(p)] - 1} 

J 47r 

The matrix element above can be evaluated and, to leading order in a Euler 
McLaurin expansion, yields p^/X^- For the phase shifts of type 6,y, associated with 
the fully interacting 3 particles, V equals V^(12) + V^(13) + ^(23) and we obtain as 
its contribution to 6f°'*^: 

ol/2 ^ 

-^Jdidff{eM-PiV{l2)+V{13)+V{23))]-l) 

The expression above, derived solely from the contribution of the 5's, diverges for 
infinite volume. However, including the terms in S, associated with the pairs 12, 13 
and 23 provides a convergent answer. The complete result for 6^°'*^ divided by b^, 
where bi — Xt, equals 

^ J dn df2 dfs {exp[-/3(y(12)+T/(13)+l/(23))] 
-exp[-/3y(12)] - exp[-(3V(13)] - ex.p[-(3V{23)] + 2} 

where 1 have integrated over R the center of mass coordinate, divided by V, and 
changed to the coordinates ri, r2andr3. The result is the classical expression with 
all the correct factors. 

Bound States 

If there are bound states, the major change in the eigenpotentials is that for some of 
these potentials, instead of going to zero at large distances (large p), there appears 
a negative 'plateau'. I.e. the eigenpotential (up to some contribution in 1/p^), 
becomes fiat and negative. This is the indication that asymptotically the physical 
system consists of a 2-body bound state and a free particle. The eigenpotential may 
also 'support' one or more 3-body bound states. 



6 



The eigenfunction expansion of the trace associated with H^, will read: 

[•OO [• r2 Q 



The q's are defined by — — 62,1, where €2,i is the binding energy of the cor- 
responding bound state. The limit Qi equals ^J^^i- The new continuum term 
represents solutions which are still oscillatory for negative energies (above that of 
the respective bound states). 

Assume, now, that we have 1 bound state, and introduce amplitudes. The 
asymptotic behaviour will be as follows. 
For E>Q. 

fAp) {kpf^Ci^i [cos 5i JKi+2{kp) - s\n5iNKt+2{kp)] 

4>\Xp) {kpy^^Ci^,,i [cos6i Jki^^Mqp) - sin^i A^i^,g+2(gp)] 
Using our procedure as before we obtain for the integral over p: 

1 r I Pmax ( sr^ 1/7 1 2 , 1/7 1 2 ^ \ 
TTdk TT q 

For E <0. 

Kip) (qpy^^icos Si JK,^+2{qp) -sin 6iNK,^+2{qp)] 



which then yields 



1 r , Pmax 
-Oi + 



TT dk Tl 

The problem is that 1 can no longer eliminate the Pmax term by subtracting the 
contribution of the free particle term; i.e. using the Pmax from T3 to cancel the Pmax 
from H-^. All is not lost however, as we saw (for example in the terms arising 
in the classical limit) that all the terms of the cluster (63) are needed to obtain a 
volume independent and convergent result. The obvious terms to examine are the 
ones associated with H2 + Ti, which also have amplitudes that correspond to (2- 
body) bound states. 1 have not been able, to date, to prove that all the coefficients 
are such that the final coefficient of pmax is zero. 

If we were ... to assume that the terms in Pmax do indeed cancel, then we can 
write the following formula for the complete trace. 

TT ^ Jo CLK 

+ ^Ee^" /%^[5,(g) -35,(?)]e-^(^^^) 
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